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ABSTRACT 

Solitary  and  periodic  internal  waves  are  shown  to  exist  in  swirling 

flow.  Incompressible/  inviscid  fluids  in  a  ri^t  cylinder  of  infinite  length 

/  . 

and  finite  radius  are  considered.  Variational  tediniques  are  used  to 
demonstrate  that  the  Euler  equations  possess  solutions  that  represent 
progressing  waves  of  permanent  form.  Moreover/  Internal  solitary  wave 
solutions  are  shown  to  arise  as  the  limiting  forms  of  internal  periodic  waves 
as  the  period  length  becomes  unbounded. 
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SIGNIFICANCE  AND  EXPLANATION 


The  Study  of  vortex  breakdown  gives  rise  to  an  interest  in  waves  in 
swirling  flow.  Our  interest  is  centered  on  the  existence  of  both  solitary  and 
periodic  internal  waves.  In  this  report  a  model  physical  problem  is  studied 
in  a  mathematically  exact  formulation.  We  restrict  our  attention  to  an 
incompressible,  inviscid  fluid  swirling  through  a  right  cylinder  of  infinite 
length  and  finite  radius.  Our  theory,  which  is  not  restricted  to  small 
amplitudes,  predicts  both  waves  of  elevation  and  depression,  depending  on  the 
angular  velocity  (swirl)  distribution  and  the  velocity  distribution  at 
infinity.  Just  as  for  the  classical  surface  solitary  waves,  these  internal 
solitary  waves  are  single-crested,  symmetric,  and  decay  exponentially  away 
from  the  crest.  Hence  they  represent  disturbances  of  essentially  finite 
extent.  Variational  techniques  and  the  theory  of  rearrangements  are  used  to 
demonstrate  these  qualitative  features.  Moreover,  we  show  that  the  solitary 
internal  wave  arises  as  a  limit  of  periodic  internal  waves  of  increasing  wave 
lengths.  . 


The  responsibility  for  the  wording  and  views  exp'^essed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


SOLITARY  AND  PERIODIC  WAVES  IN  SWIRLING  FLOW 


Scott  A.  Markel 


1.  Introduction 

The  Interest  In  waves  of  finite  amplitude  and  permanent  form  in  swirling  flow  arises 
from  the  study  of  vortex  breakdoim.  This  is  the  rapid  change  in  structure  which  can  occur 
in  swirling  flow.  An  example  of  this  is  vortex  breakdown  above  a  triangular  wing.  The 
reader  is  referred  to  the  book  by  Van  Dyke  (1982)  which  contains  excellent  photographs  of 
vortex  breakdoim. 

We  restrict  our  attention  to  an  incompressible,  Invlscld  fluid  swirling  through  a 
right  cylinder  of  infinite  length  and  finite  radius.  The  waves  studied  in  this  paper  are 
the  analogues  in  a  swirling  flow  of  the  internal  solitary  and  periodic  gravity  waves 
discussed  in  Bona,  Bose,  Turner  (1983).  For  this  reason  we  will  call  them  solitary  and 
periodic  waves.  The  variational  approach  used  here  closely  follows  the  approach  used  in 
Bona,  Bose,  Turner  (1983).  Their  introduction  contains  a  survey  of  the  literature  on 
Internal  and  solitary  waves.  The  papers  of  Benjamin  (1962)  and  Pritchard  (1970)  discuss 
the  theory  of  vortex  breakdown.  Pritchard  also  descrl)>es  various  experiments  in  rotating 
flows  in  which  solitary  waves  were  observed. 

In  section  2  the  Idealized  physical  model  is  described  and  the  governing  equations 
are  derived.  The  main  item  of  interest  in  this  paper  is  the  apparent  difficulty  arising 
from  the  singularity  of  the  governing  equations  along  the  axis  of  the  underlying 


cylindrical  domain.  The  existence  theory  can  be  readily  adapted  to  cover  the  singular 


case.  The  existence  results  for  Internal  periodic  waves  are  given  in  section  3.  Both  the 
case  where  the  speed  of  propagation  Is  specified  and  the  case  where  the  wave  energy  Is 
specified  are  considered.  Since  the  governing  equations  are  singular  along  the  axis  of 
the  cylindrical  domain,  standard  elliptic  theory  is  not  applicable  there.  By  restating 
the  problem  in  a  higher  dimensional  space,  the  desired  regularity  of  solutions  is 
achieved.  This  is  proved  in  section  4  and  uses  an  idea  due  to  Ni  (1980).  Section  5 
begins  with  a  discussion  of  a  priori  bounds  satisfied  by  the  periodic  solutions.  These 
bounds  are  independent  of  the  period  length  and  are  found  in  section  4  of  Bona,  Bose, 
Turner  (1983).  A  method  due  to  Amiclc  (1984)  is  used  to  obtain  the  exponential  decay  of 
periodic  waves  from  crest  to  trough.  The  final  result  of  section  5  is  the  existence  of 
internal  solitary  waves.  They  are  shown  to  arise  as  the  limiting  forms  of  internal 
periodic  waves  as  the  period  length  becomes  unbounded.  A  specific  exanple  of  swirling 
flow  is  discussed  in  section  6  and  is  shovm  to  yield  governing  equations  which  fall  within 
the  limits  of  our  theory. 

The  author  wishes  to  thank  Prof.  Robert  E.L.  Turner  for  his  many  helpful  suggestions. 


2.  The  Governing  Equations 

To  investigate  the  existence  of  internal  waves  in  swirling  flow,  an  idealized 
physical  model  is  considered.  Attention  is  restricted  to  an  incompressible,  inviscid 
fluid  swirling  through  a  right  cylinder  of  infinite  length  and  finite  radius.  Two- 
dimensional  flows  will  be  our  domain  of  interest  by  assuming  the  flows  to  be 
axisymmetric.  We  will  use  cylindrical  coordinates  with  the  z-axis  as  the  axis  of  the 
cylinde’-  and  r,  the  radial  distance  from  the  center  of  the  cylinder.  The  radius  of  the 
cylinder  is  ta)cen  to  be  a,  and  thus  the  boundary  is  given  by 
( (a,9,z)  :  0  f  [0,27tl  ,  z  «  r}. 
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where  '  denotes  differentiation  with  respect  to  i|i.  Equation  (2.5)  is  the  fundamental 

equation  describing  our  stodel  problem.  The  corresponding  kinematic  conditions  are 

i|i(0,z)  >  ftO)  -  0  and  i|>(a,B)  -  y(a),  ^  (2.6) 

for  z  c  It  and  the  asymptotic  condition  Is 

4i(r,z)  "Kr)  as  jz]  ♦  “  (2.7) 

for  r  e  [O.a,].  i)i  represents  a  flow  connected  to  the  primary  flow  at  Infinity  and  for 

which  I  and  H  are  constant  on  stream  surfaces.  H  and  I  are  not  immediately  known 

as  functions  of  i|/,  however  In  principle  they  may  be  determined  from  the  primary  flow. 

—  2 

For  exaiq>le.  If  W(r)  =  d,  a  constant,  and  c  •  d  -  c,  then  by  (2.2)  y(r)  •  cr  /2 

2^  1/2 

which  may  be  Inverted.  The  Inverse  Is  RCf)  •  (~)  *  'circulation"  I  Is  then 

expressed  as  a  function  of  the  stream  function  value  V  by  1(7)  -  i(r(7)).  Since  I  Is 

constant  on  stream  surfaces.  Its  value  at  a  point  flow  corresponding 

to  i|i.  Is  determined  by  tracing  the  stream  surface  with  value  infinity 

where  I  can  be  determined  by  I(V)  •  i(r(7)).  H  Is  determined  from  H  -  —  t  — |q|  and 

2 

the  equation  of  radial  equilibrium  In  a  steady  cylindrical  flow  :  ~  pv  /r.  This  Is 

obtained  from  (2.3).  See  the  exaiq>le  In  section  6  for  details  on  how  to  determine  1(7) 

and  H(7)  for  a  specific  problem. 

Define  the  perturbed  stream  function  (p(r,z)  by 

i(i(r,z)  -  7(r)  +  y9(r,z),  (2.8) 

where  y  Is  a  normalizing  constant.  In  the  primary  flow  we  have 

7„  -  -  7  -  r^H’(7)  -  I' (7)  . 

rr  r  r 

Hence,  the  form  of  H'  is 

H'('f)--7^ -  r7^^(7)  -  7  (7)  +  I'(7)l  . 

R^(7) 

By  following  stream  surfaces  to  Infinity  and  remembering  that  H  Is  constant  on  stream 


surfaces,  we  find  that 


■  --  ^  ^  ^ 


•  ■r 


jrjrj^j»i»pj  iji  jijjiuulu 


1f(9  ”  “  <P  *9  )  “  I'(T  ♦  Y?)  [ 

I'-Tpp  r  fj,  Tjj  I 


R  (T  +  Y9) 


R^{t) 


R  (T  +  Y9) 


%(»  +  Y») 


-  t’-rr'”  -  %r]-  «• 

For  further  details  on  the  derivation  of  (2.9)  see  Benjamin  (1971),  where  an  anal090us 
problem  in  stratified  flow  la  considered. 

Equation  (2.9)  can  be  written  as  the  nonlinear  eigenvalue  problem 

-Ls  +  h(r,f)  “  Xf(r,e) 


where 


r  2r 


3s 


(2.10) 


The  constant  X  is  an  eigenvalue  perimeter  proportional  to  c~^ ,  where  c  is  a  velocity 

scale  for  the  primary  flow  in  the  travelling  coordinates  introduced  earlleri  eg. 

c  >  sup  H(r).  The  supplementary  conditions  are 
0<r<a 


s(0,t)  -  v(a,B)  -  0,  z  €  R, 
v(r,z)  >0  as  ts|  r  €  [0,a]. 


(2.11) 


In  this  paper  it  will  be  shown  that,  under  suitable  conditions  on  f  and  h,  the 
boundary  value  problem  (2.10,  2.11)  has  solitary  wave  solutions.  By  a  'solitary  wave" 
solution  we  mean  a  solution  9  which  is  even  in  z,  monotone  for  z  >  0,  and  rapidly 
convergent  to  zero  as  |z|  In  the  course  of  obtaining  solitary  wave  solutions  it 

will  be  shown  that  there  are  solutions  9  of  (2.10)  which  vanish  at  r  •  0  and  r  *  a, 
which  are  even,  periodic  functions  of  z  with  period  2R,  and  which  are  monotone  for 
z  c  [0,k].  These  will  be  called  periodic  waves. 

3.  Periodic  Waves 

In  this  section  we  prove  the  existence  of  periodic  solutions  (1,9)  of  the  boundary 
value  problem 
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-LfCr.z)  +  h(r,9(r,z))  -  Xf(r,9(r,z))  in  n, 
^(0,z)  -  7(a,z)  «  0, 


(P)  (3.1) 


where  L  is  given  by  (2.10)  and  fi  »  ((r,z)  s  r  c  (0,a),  z  £  r}. 
He  make  the  following  hypotheses  on  f  and  h. 


Hypothesis  (H); 


(HI):  The  function  f  has  the  form 


f(r,t) 


t  +  f^(r,t)  r  €  (0,aj,  t  >  0 


-  f(r,  -  t) 


r  €  I0,a)  ,  t  <  0 


with  ^o'^l  Hdlder  continuous  on  bounded  r  sets  and  f^  Lipschitz  continuous  in  t  on 
bounded  sets.  He  further  assume  that  fg  >  0  for  r  £  (0,a),  f^(r,t)  =  o(t),  uniformly 
for  r  £  (0,a],  as  t  «  Ot  and  that  there  exist  constants  m,n  >  1  and  a.d  >  0  such 


that  for  t  >  0 


at“  <  f,(r,t)  <  d(1  +  t”). 


(HID:  The  function  h  has  the  form 


h(r,t)  = 


t  hg(r)  +  h^(r,t)  r  £  [0,a) ,  t  >  0 

-  h(r,  -  t)  r  £  [O.a],  t  <  0 


with  hg,hf  Hdlder  continuous  on  bounded  r  sets  and  h^  Lipschitz  continuous  in  t 

d2  ,  j 

on  bounded  sets.  Let  eg  be  the  lowest  eigenvalue  of  -  — —  ~  in  [0,a]  with 

dr 

Dlrichlet  data.  A  lower  bound  for  eg  is  given  in  lemma  3.4.  For  h  we  further  assume 

that  is  continuous,  —  i  e  where  e  >  “  e  ,  and  that  there  exist  constants  o  »  0 

dt  0 


and  c'  such  that  for  t  >  0, 


c't  '  (r,t)  <  ot  , 


|h  (  r,  t )  I  <  d( 1  +  t") , 


where  Bin,  end  d  are  the  constants  In  (HI)*  The  first  Inequality  Implies  h^(r,t)  • 
o(t),  uniformly  for  r  e  [OiS],  as  t  >  0* 

Finally  ire  define  the  functions 


t  t 

P(r,t)  -  2/  f(r,s)d8,  F  (r.t)  -  2/  f  (r,s)ds 
0  0 

t  t 

H(r,t)  -  2/  h(r,s)ds,  H  (r.t)  «  2/  h  (r,8)ds. 
0  0  ’ 


(3.2) 


(Hill):  The  last  hypothesis  Is, 

For  each  X  >  a/ a,  where  a  and  0  appear  In  (HI)  and  (HID,  there  Is  a 
e  «  (0,1)  such  that 

XF^(r,t)  -  H^(r,t)  <  e{xf^(r.t)  -  h^(r.t))t. 

We  now  want  to  find  solutions  of  (P)  which  are  periodic  In  e  with  period  21c. 


Let 


2  "  l(t.a)  8  r  €  (0,a),  z  e  (-k.k)}.  (3.3) 

The  subscript  2  here  is  to  distinguish  the  set  from  higher  dimensional  analogues  to  be 
Introduced  later.  The  problem  (P)  can  be  formulated  In  two  ways.  The  first  is  a 
constrained  problem  (PC): 

solve  ( P ) 


12  2 
subject  to  /  ~  {1^9!  *  H(r,9)}  drdr  «  R 

«k.2 


where  R  >  0  Is  a  given  constant. 

The  second  is  a  free  problem  (PF): 

solve  (P) 

where  X  is  a  given  constant. 


(PC)  (3.4) 


(PF)  (3.5) 


(PC)  corresponds  to  specifying  the  "energy"  of  a  wave,  while  (PF)  corresponds  to 
specifying  Its  velocity. 

The  analysis  of  both  problems  Is  based  on  variational  techniques  and  the  theory  of 


rearrangements  of  functions.  In  each  problem  a  critical  point  9  of  a  suitable 


functional  will  be  found  and  shown  to  be  a  solution  of  a  weak  formulation  of  the 
problem.  In  addition,  9  can  be  taken  to  be  even  in  z,  nonnegative,  and  ^tonincreasing 
on  z  €  [0,k)  for  each  r.  The  following  definition  is  used. 


Definition  3.1;  Let  »  “  »(r,z)  be  continuous  on  2  *  with 

period  2k.  For  each  r  let  ii(T,T<r)  denote  the  Lebesque  measure  of  the  set 
{z  :  9(r,z)  >  y).  A  function  9  which  is  even  in  z,  nonnegative,  nonincreasing  on 
z  €  (0,k]  for  each  r  and  satisfies  11(9, Y**^)  “  w(|9|,Y»t)  is  called  the 
symmetrization  of  9.  If  9  “  9  we  call  it  symmetrized. 


The  analysis  will  be  carried  out  in  the  Hilbert  space  =  H^((0,a)  x  r]  defined 

as  follows:  let  C*  denote  the  C°°  functions  which  have  support  where  r  c  (0,a)  and 
k 

which  are  periodic  in  z  with  period  2k.  Define 

lul^  =  /  ^  {1^1^  drdz.  (3.6) 

”k,2 

The  Poincare  inequality  (see  lemma  3.4) 

f  ~|?u|^drdz  >  e„  /  “  u^  drdz  (3.7) 

\,2  "  «k,2  " 


together  with  ®  ^  ~  ®0'  *  consequence  of  (HII),  show  that  *  is  a  norm. 

«0 

Let  H,  be  the  completion  of  C  in  this  norm;  the  space  H.  is  a  Hilbert  space  with 

K  K  K 


inner  p'-oduct 


(u,v)k 


/  ^  (Vu«Vv  +  hp(r)uv}  drdz. 

\,2 


(3.8) 


If  I  is  a  continuous  linear  functional  on  Hj^,  its  value  at  u  is  denoted  by  <f,u>. 

The  use  of  synmietrized  functions  is  not  needed  in  showing  the  existence  of  periodic 
solutions,  however  it  is  very  important  in  getting  the  estimates  in  section  5  which  are 
used  to  show  that  a  solitary  wave  can  be  obtained  as  the  limit  of  periodic  waves  with 
increasing  period.  For  the  construction  of  9,  using  piecewise  linear  functions  and  an 


approximation  process,  and  the  properties  of  9 


the  reader  is  referred  to  the  appendix 


of  Bona,  Bose,  Turner  (1983).  A  version  of  the  following  lemma  can  be  found  In  Polya, 
Szego  (1951). 


Lemma  3.2;  Suppose  G(r,u)  Is  even  In  u  and  continuous  for  (r,u)  e  (0,a]  x  r.  Then, 

for  u  piecewise  linear  and  u^  0  a.e., 

/  G(r,u)drdz  •  /  G(r,u)drdz.  (3.91 

\,2  \,2 

Lemma  3.3;  Suppose  p(r)  Is  positive  and  measurable  on  [0,a)  and  u  Is  piecewise 

linear  and  21c  periodic  In  z  with  u^  0  a.e.  Then 

/  p(r ) I Vuj ^drdz  >  /  p(r) I Vul^drdz.  (3.10) 

\,2  \.2 


Proof ;  The  case  p  z  1  Is  found  In  Polya,  Szego  (1951).  By  continuity  we  may  assume 

p  «  C^([0,a])  and  u  c  C^.  Let  “  p  with  q  (r)  >  0.  Consider 

2  2  2  2  2 
|7(qu)|  •  q  |7ul  ♦  q^u  2qq^uu^. 


/  qyq(u^) jdrd*  “  -/  (q^^qu*  ♦  q*u*)drdz. 


2  2 
/  p|Vu|  drdz  -  /  |q7u|  drdz 


/  |7(qu?|  drdz  -  /  (q^u^  +  2q  qu  u)drdz 

«lc,2  «k,2  " 


2  2 

-  /  l7(qu)|  drdz  /  **rr'*'*  drdz. 


Since  the  symmetrlzatlon  Is  In  z  and  since  q  is  a  positive  function  of  r,  qu  “  qu. 

By  lemma  3.2,  with  G(r,u)  -  qyj(r)q(r)u^(r,z), 

2  *'2 
/  '^rr*^  drdz  «  f  q^^qu  drdz. 
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Finally 


/  pI^I  drdz  «  /  Ivtqu)!  drdz  +  /  q  qu^drdz 


^,2 


”k,3 


V2 


>  /  |7(qa)|  drdz  +  /  q  qu  drdz 


«  /  |?(qu)|  drdz  +  /  q  qu  drdz 

«k.2  «)c.2  "" 


/  p|7u|  drdz* 

”lc,2 


Lemma  3.4:  For  u  c  the  constant  eg  in  the  Poincare  inequality  (3.7)  satisfies 

4ii^+3 
e.  >  - :r“* 

°  4a2 


Proof ;  By  continuity  we  assume  u  £  Cj^.  Using  the  calculation  in  the  proof  of  lemma  3.3 


with  q  =  1//r, 


f  ^  |7u|^drdz  »  /  |7(^)l^drdz  +  /  drdz 

'  a  -  4r^ 


^,2 


2  2  ,2 

>  ^  ^  drdz  +  — r  /  ^  drdz 

a  Cl.  .  4a 


r  ^  u^  drdz. 
«K,2" 


Along  with  the  spaces  C^,  L**(n),  and  W^'^(n)  will  be  used  for  various 

domains.  The  reader  Is  referred  to  Adams  (1975)  or  Gilbarg,  Trudinger  (1983).  Here  and 
In  what  follows  C  denotes  a  positive  constant  independent  of  k.  The  constant  changes 
from  one  inequality  to  another. 


Lcawa  3«St  The  following  relation  holds « 


p  e  I2,«). 


Proof;  From  the  Poincare  inequality  and  Sobolev  imbedding  theorem, 
lul  _  “  (/  lupdrdz)^^*^  2  <  p  <  •• 


”k,2)  \,J 


<  C  (/  {1^1^  +  u^}drd*)'^^ 

-  C  lul 

«  '  ‘V2> 

<  C  /a(/  —  {|Ptt|^+  u^}drd2)^^^ 

<  C  (/  ^  {|Pu|^+  h  u^ldrdm)’''* 

”lc,2 

■  Clul.  . 


A  further  consequence  of  the  Poincare  inequality  is,  for  u  e 


lul^  <  C  /  |Wi|^drds. 

”)t,2 


For  convenience  we  will  often  use 


/u  “  /  u(r,z)drdz. 

”lt,2 


If  u  €  Lu//r  «  2*  •"**  ''  ^  "k' 


(u,v)|^  "  ^  (Tu'Vv  +  hp(r)uv) 


(3.11) 


(3.12) 


(3.13) 


where  the  boundary  terms  cancel  on  a  <■  dt  k  due  to  periodicity.  Thus  %fe  define  a  weak 


periodic  solution  of  (P)  as  a  pair  (X/f)/  e  e  satisfying 

(».v)k  +  /  ^  h^(r,,)w  -  Xf  j  f(r,9)v  (3.14 

for  all  V  e  H^.  The  Integrals  In  (3.14)  exist  fay  virtue  of  (HI),  (HID,  lemma  3.5,  and 
lemma  3.6. 


Lemma  3.6;  For  »  «  P  >  2, 


f  ^lul**  <  C  lul^. 


(3.15 


Proof;  Let  q  -  1/p,  then 


/  il-l”  - 1  -  fcl 


From  the  proof  of  leana  3.5, 


<  c 


1/2 


Since  p  >  2,  q  <  1/2  and  r"***  <  C  r”'. 

Considering  now  the  gradient  term  and  Integrating  by  parts  as  in  the  proofs  of  lemmas  3.3 
and  3.4, 


<e  <,%.  ;  J 


2q+2 


Hence 


/  /  {^\f}  4C  .U. 


2  2 


and  the  result  follows. 


First  we  consider  the  problem  (PC)«  For  u  c  introduce 
A(u)  “  /  ^  {|7u|^+  il(r,u)}  'I 

J(u)  “  /  ^  F(r,u)  I 

where  F  and  H  are  defined  in  (3>2)  and  for  R  >  0,  let 

S(R)  "  (u  «  :  A(u)  -  R*}.  ^ 


C3.1 


The  aasui^tlons  on  f  and  h,  along  with  lemma  3>6,  guarantee  that  the  functionals  A 
and  J  are  defined  on  One  calls  s  a  critical  point  of  J  on  S(R)  if  the 

derivative  of  J  la  saro  in  directions  tangent  to  6{R),  le.  J'(s)  is  parallel  to 
A'<s).  A  consequence  of  this  la  (3.14). 


Theorem  3.7;  Let  f  and  h  satisfy  (H>.  Tfien  for  each  k  >  0  (PC)  has  a  solution 

(Xj^.e^).  #1^  t  ’Sc'^ 

(1)  J(e.  )  ■  sup  J(u), 

u  (  8(R) 

(2)  >  0  and  fjj  >  0  in  ^  , 

(3)  (cf.  definition  3.1). 

Proof !  Apart  from  regularity,  the  proof  follows  the  lines  of  theorem  3.2  in  Bona.  Bose. 
Turner  (1983).  Due  to  the  1/r  singularity  in  the  operator  L  (2.10).  elliptic 

theory  from  Agmon,  oouglis.  Rlrenberg  (1959)  is  applicable  only  to  regions 
n'  ~  {(r.z)  c  2  I  r  >  <  >  O}.  8ince  y  is  a  weak  solution  of  an  elliptic  equation. 

2,p 

e  €  H  (O'),  p  c  [2.»).  By  the  Bobolev  imbedding  theorem  and  the  Schauder  theory,  see 
Gilbarg,  Trudlnger  (1983),  e  €C*'“<n*). 

Remark;  In  section  4  it  is  shown  that  g  «  j),  ie.  that  e  regular  at  r  « 


Now  we  turn  to  (PF).  Let  u  be  the  smallest  eigenvalue  of  the  linear  problem 


(3.17) 


■Lu  +  hgCrlu  -  Xfjjlrlu, 


u  € 

Using  separation  of  variables  we  see  that  p  has  a  corresponding  eigenfunction  v(r) 
associated  with  the  lowest  eigenvalue  of 

-''rr  *  r  V  *  Xf^CDv. 

(3.18) 


v(0)  “  u(a)  “  0. 

3 

In  fact  the  equation  is  of  limit  point  type  at  r  -  0  with  respect  to  L  with  the 
weight  1/r,  so  the  condition  v(0)  «  0  is  superfluous.  We  can  also  use  the  change  of 

3 

variables  y  -  r  /2  to  get 


5'g(y>v  Xfplylv 

■’'yy  ^  27“' 

3 

where  v{y)  v(r  /2)  »  v(r).  The  existence  and  positivity  of  v  follow  from  arguments 
like  those  used  in  the  proof  of  theorem  3.7.  We  normalize  v  such  that 

/  i  fv^  +  h  v^}dr  •  1.  (3.19) 

0 

For  X  fixed  define  a  functional  M  on  by 

M(u)  -  lul^  /  r  (H,(Jf*u)  -  XF(r,u)}.  (3.20) 

IM'(u)l  denotes  the  norm  of  the  derivative  M' (u)  as  a  functional  on  Hj^.  A  critical 
point  ®  of  M,  ie.  <M'(?>),v>  -  0  for  all  v  e  is  a  weak  solution  of  (PF)  l>y 

(3.14).  Now  we  show  (PF)  has  a  symmetrized  solution  by  first  stating  a  technical  result 
from  Bona,  Bose,  Turner  (1983). 


Proposition  3.8;  Fix  k  >  0  and  let  H  be  a  continuously  differentiable  functional  on 
satisfying: 


(1) 

There  are 

constants 

s,  s  >  0 

such  that 

N(u)  >  s  for 

lull,  =  s. 
k 

(2) 

M(0)  =  0 

and  there 

exists  a 

function  w 

with  Iwl,  >  s 
k 

such  that 

M(w)  <  0. 

(3) 

For  each 

B  >  0  if 

00 

is  a  sequence  satisifylng 

B  <  M(u^) 

<  1/B,  for 

1, 

and  HM ' (u 

)  n  0  as 

1  ®, 

then  there 

is  a  subsequence 

(“i  Cl 

of 

m 


eo 

converging  strongly  in  as  m  +  “. 
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r  -  {Y  e  C([0,1],Hj^)  .  YlO)  -  0,  Yd)  -  W) 


and  define 


b  -  b„(M) 


u  «  Y( 10,1) > 


b  -  inf  bY* 

Y  «  r 

Suppoae  Y  •  n  c  ia  a  sequence  of  paths  In  r  auch  that  b  «  b  'f  1/n.  Then  there 

Is  a  subsequence  {n.}.  ,  of  S  and  functions  u  e  Y  (t0,11)  such  that  u_ 

j  n^  n^  "j 

converqes  stronqly  to  u  as  j  M(u)  -  b,  and  M'(u)  •  0. 


Theorem  3.9i  Let  f  and  h,  satisfy  (H)  with  a  and  o  as  defined.  Let  u  )>•  the 
lowest  eigenvalue  of  (3.17).  if  X  e  (a/a,u)f  then  for  each  k  >  0  (PF)  has  a  solution 


’k*  '  <«k.2» 

(1)  M(e  )  -  inf  b  (M), 

Y  <  r  ^ 

(2)  e^>0  in 

d) 


Proof »  The  proof  parallels  that  of  theorem  3.4  in  Bona,  Bose,  Turner  (1983).  One 
constructs  a  sequence  of  minimising  paths  consisting  of  symmetrised  functions. 
Proposition  3.8  guarantees  that  the  critical  point  will  be  symmetrised. 


Remark!  Again  we  postpone  the  regularity  of  9  at  r  ~  0  until  section  4. 


4.  Regularity  At  r  ~  0 

In  this  section  we  obtain  regularity  up  to  r  ~  0  where  the  equation 


-(9 _ ■  Z  7 _ )  h(t»e)  "  Xf(r,9) 


is  singular  and  standard  elliptic  theory  breaks  down.  To  extend  the  regularity  of  f 


from  to  we  use  an  approach  due  to  Hi  (1980),  which  Involves 


restating  the  problem  in  H 


Set  9  «  r  g«  We  will  fix  k  for  the  entire  section  and  suppress  it.  We  then  have 


Lo  "9  •  “  9  +9 

V  y  ▼jp. 


r^grr  +  3rg^  + 


Henc«  (4.1)  becomes 


'■’rr  r  ’t  ’zs-'  _2  2 


(4.2 


with  boundary  conditions 

g(a,z)  •  0,  g(r,k)  -  g(r,-)c). 

The  operator  which  has  replaced  L  is  the  LapXaclan  in  cylindrical  coordinates  in  a 

applied  to  functions  with  cylindrical  symmetry.  Since  the  nonlinear  terns  f  and  h 

have  polynomial  growth  in  (r  g)<  division  by  r^  is  not  a  problem. 

We  define  g  in  to  be  axially  symmetric  with  respect  to  the  axis  r  •>  0  (z- 

axis)  such  that  g|  „  **  V*  (V,  e  domain  in  generated  by  rotating 

'‘>,2 

around  the  axis  r  ~  0.  Note  that  r  •  0  is  not  included  in  ^5* 


”h,2 


Lemma  4.1;  g  is  a  classical  solution  of  (4.2)  in  the  open  set  g. 


Proof;  This  follows  Immediately  from  knowing  ?  -  r  gjjj^  ^ 
(3.1)  for  r  >  0. 


is  a  classical  solution  of 


We  will  use  j)  and  g)  to  denote  the  and  w^'P  spaces  in  R®.  We 

omit  the  angular  factor  of  the  Jacobian,  r^  aln^e^  sine^  drde^  de^de^dz,  in  the 


integrals  for  n 


Lftinma  4.2:  For  f  c  and 

•Sflfc  ■  /  ll’9l^  •*•  h  g^}r^drd*  -  l»l^. 

Vs 


Proof:  IV(r^g)|^ 

»  r^  I  7gl  ^  t  4r 

3  ,22 

99j.  +  4r  g 

and 

/ 

9  J  2 

drd*  “  ”  /  rg  drds. 

in 

r 

so 

/  7I 

7(r^g)l^drds  - 

/  r^|7g|^drds. 

Therefore, 

"k.2 

Vs 

•»lk  -  /  J  {lv,|^  +  h^(r),^}  drd* 

V2 

-/  (|7g|V  (r^g)^}  drd* 

"k,5 

■  /  {|Vg(^  +  h-g^}  rVdr 

"k.5  , 


Theorem  4.3i  g  Is  a  weak  solution  of 


Xf(r,r^q)  -  h(r,r^g) 


with  g(a,z)  “  0  for  *  e  (-k,k)  and  g  periodic  In  z  with  period  2k.  ®)t,5  “ 

{(O.z)  ;  z  e  (-k.k)}  U  0^  5. 

Proof!  The  reader  la  referred  to  N1  (1980),  where  It  Is  shown  dg  does  not  produce  a 


distribution  supported  at  r  •  0. 


Iieima  4.4;  Suppose  u  -  u(r,z)  t  W  '  (B  )  with  u  (a,*)  ■  0  for 
~  #  5 

3q 

z  €  (-k.k),  q  >  1,  and  6  >  j*  Then 


(/  |u|**r®drdz)^'^^  <  c(/  |  Vu|^r^drdz)^^^. 


Proof ;  By  continuity  vre  assume  u  «  c”(B  ).  The  Sobolev  inequality  in  is 

k  1 5 


Consider 


(/  |ur®^‘*|'*drdz]’/‘*  <  C(/  |7(ur®^‘*)|^drdz)'^^ 


,  6_  ^2  ^(2B/q>  -  2^ 

q 

,  8,2  (26/q)  -  1  .  .  6  26x  f  2  (28/q)  -  2  .  ^ 

\,2  \,2 

r  |7(ur*’  '‘*)f^drdz  -  /  |  7u|  ^^^^'^drdz  +  ^  d"  ^  drdz 


<  /  |7u|^^®^Vdz 

V2 


,  (2B/q)  -  3  f  In  |2  3.  _ 

<  sup  r  j  I 7u I  r  drdz 


<  C  /  |7ul^r®drdz. 


The  result  follows  immediately. 


Now  we  state  the  main  res-ilt  of  this  section. 


Theorem  4.5;  g  •  or  equivalently  f/r^  t  ^) , 


1.'  ■ 


*^e  “^e  “^e 


•  .'•'.••V-V-V- 


J 

...a 
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••  v 


Proof t  By  (H),  |xf(r,t)  -  h(r,t)|  <  C(|t|  |t|''),  and  so  the  right  hand  side  of  (4.3) 

la  bounded I 


Xflr,r  a).  -  htr.r  al|  «  <.(1,1  ♦  r^’‘"*|,r). 


Standard  elliptic  theory  la  applicable  to  (4.3)  since  we  are  now  dealing  with  the 

1  2 

Laplaclan.  Xieana  4.2  shows  g  c  H  '  (B  .).  For  n  <  7/3,  a  bootstrapping  arguswnt  and 


the  crude  estlnate  C(|g|  +  r*""^|g|")  <  C(|gj  +  |gl")  give  us  g  e  w^'^(B.  _)  for  all 

K/5 

p  >  1.  Invoicing  leana  4.4  with  p  -  5/2,  q  •  np,  and  6  ~  (2n-2)p  *■  3, 

,  ,  2n-2|  ,n,p  3.  .  , 

J  (r  |gl  )  r  drds  <  <• 

®k,5 

If  (2n-2)  '2  +  3  *  2^  ’  ^  "  *  B/i,  Hence,  for  n  >  8/5,  the  right  hand  side  of 


(4.3)  Is  In  l‘^^(B]j,s)*  Thus,  by  elliptic  theory,  get)' 

p 

follows  free  the  Sobolev  Inbeddlng  theorea  that  g  c  L  (B  .)  for  p  e  [I,**).  Since,  by 

Itp  5 

2sP 

elliptic  theory^  g  €  W  (B.  «)  for  p  €  a  further  application  of  the  Sobolev 

K/5 


2,5/2 


(B  .)  for  all  n  >  1.  It 

JC  y  5 


Inbeddlng  theorea  gives  g  e 

,  c“(B,  ,). 


K,5' 


(B.  ...  By  our  aasuaptlona  on  f  and  h, 

K  p  d  ) 

The  Schauder  theory  gives  g  e  C^'“(B  ).  Our  result 

JC  p  o 


follows  froa  noting  that  B 


k,5 


‘’lt,5* 


5.  Solitary  Waves 

As  In  Bona,  Bose,  Turner  (1983)  It  can  be  shown  that  <  W  end  bounded 

Independent  of  k  for  both  (PF)  and  (PC).  The  additional  assuaptlons  of  a  <  5  and 

2a/a  <  u  ere  needed  here,  where  a, o,  and  a  appear  In  (H).  Corollary  4.5  In  Bona, 

Bose,  Turner  (1983)  gives  a  bound  on  Ig  1  .  Independent  of  k.  Since  our  solutions  of 

k  ^2,0. 

(PF)  and  (PC)  ere  syaaetrlsed  In  z,  for  k  sufficiently  large  we  can  obtain  the 
following  crude  decay  estimate 


g^{r,«) 


'='%')c 


(1 


1 1/3 


(5.1) 


for  |s|  <  k.  See  lemma  4.9  of  Bona,  Bose,  Turner  (1983)  for  details. 
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Osing  this  dscay  rats  ws  can  obtain  exponential  decay  rates  for  9^  and  The 

approach  m  use  Is  due  to  Udch  (1984).  He  used  It  to  get  exponential  decay  for  solitary 
waves,  but  It  Is  easily  adapted  to  periodic  waves. 

Leana  5.1;  Let  be  a  solution  of  (P)  from  either  theorem' 3.7  or 

theorem  3.9.  Suppose  that  for  some  q  >  1, 

|)i^f,(r,t)  -  h^(r,t)|  <  Cpt^ 

for  |t|  <  1  and  r  e  [0,a].  Then  for  k  >  k^,  |z|  <  k, 

9j^(r,2)  <  Ce  (5.2 

and 

|V9  (r,z)|  <  C'e"®l'l  (5.3 

K  V2 


where  C  and  C'  depend  on  a,q,6,  and  Ig.  I 


«k'c2,<o 


e  t  e. 


,)0 


Proofs 

Q(r,g) 

and 


As  noted  above,  <  li. 
Xf^(r,r^g)  -  h^(r,r*g) 


We  now  suppress  the  subscript  k.  Let 
2 

Again  we  use  9  -  r  g|  where  g  satisfies  (4.2), 

'"k,2 


lQ(r,g)|  <  C^g**  If  r^g  <  1  . 

For  r  €  (0,a],  z  £  [-k,k] ,  (5.1)  gives  g(r,z)  <  - z-rr.  For  k.  large  enou^ 

2  I'l 

r  g  <  1  on  [0,a]  x  [kQ,k].  We  can  rewrite  (4.2)  as 

-7-(r^Vg)  +  h^gr^  -  ^^0^^  *  • 

Multiplying  by  g  and  Integrating  over  (0,a)  x  (z,k)  gives,  after  an  Integration  by 
parts, 

ka  ka 

/  /  {l’9l  +  h  g^fr^drdc  -  /  /  Xf  g^r^drd? 

z  0  z  0 


(5.4 


\i  Is  the  smallest  eigenvalue  of  the  linear  problem  (3-17)  and  g^tr^k)  <  0  since 


g  *  g,  so 

ka*-*  a  ka  * 

(■'-  ”)  /  /  {l’9r  +  h  9  }r’drdc  <  -  /  g(r,«)g  (r,*)r  dr  +  /  /  Q(r,g)gr  drdi;. 
**  *  0  0  «  0 

rroa  the  proof  of  lenoe  4.2  and  our  Poincare  inequality  (3.7) 

/  /|Vg|*r^drdc  >  e^j  /  /  g^r^drdc- 

Using  this  and  the  hypothesis  on  1q|  ,  (5.5)  becomes,  for  z  e  (Ic  ,k). 


/  /  <*o  hQ>g  r  drdc  <  /  /  g**  r^drd?  -  /  g(r,*)g^(r,*)r  dr. 

**  z  0  *  0  0 

From  (HID  h„(r)  >  e  >  -  e„,  so  let  C,  -  (e„  +  e)(l  -  — ).  Then 

0  0  3  0  W 

k  a  —  —  a  ^  ci*1  k  a  2  3 

C  /  /  g  r  drdt  +  /  g(r,«)g^(r,s)r  dr  <  max  |g(r,E  )!'*"//  ,"r"drdc. 


for  z  €  (k^.k).  Choose  kQ  such  that 

C  q-1  C 

C,  max  Ig(r,«)|‘*  <  C,( — . 


Now  for  z  e  (k.,k)  we  have 
0 


C-  ^*31  3 

—  J  J  g  r  drdc  +  /  g(r,z)g  (r,z)r  dr  <  0 


As  in  Amick  (1984),  it  follows  that 

’'“23  ' 

/  /  g  r  drdc  <  C  e 
z  0  4 

for  z  c  (k.,k).  Since  Igl  „  is  bounded,  we  obtain 

0  ^,2,(1) 

g(r,z)  <  C^e 

on  (0,a)  X  [-k,k),  k  >  k„,  where  8  “  /C,/2. 


Elliptic  estimates  (eg.  iP  estimates  and  the  Sobolev  inequalities)  give  local  l” 
bounds  for  7g  in  terms  of  local  L  bounds  for  g.  It  follows  that 

I Vg(r,z) I  <  C^e  . 

We  now  show  the  existence  of  solitary  waves  for  (PF)  and  (PC).  They  arise  as  the 
limit  of  periodic  waves  as  the  period  k  +  «>. 

Theorem  5.2:  Assume  f  and  h  satisfy  (H)  and  let  a, o#  and  m  be  the  constants 
defined  there.  Let  g  be  the  lowest  eigenvalue  of  the  linear  problem  (3.17).  Suppose 
also  that  m  <  5,  2o/a  <  g,  and  |Xf^(r,t)  -  h^(r,t)|  <  Ct**,  for  some  q  >  1/  on 
[O.a]  X  [-1,1].  Then 


(PC):  Let  R  >  0  be  given.  For  each  k  >  0  let  *  solution  of  (PC)  from 

2 

theorem  3.7  with  •  Then  there  exists  an  increasing  sequence  of  half  periods 

k(j)  j  *  Z*,  and  a  solution  (X,(p)  of  (P)  satisfying 


(1) 

CP  >  0 

and 

^  «  <p 

on  Cl, 

(2) 

.1  («) 

=  R^ 

(3) 

4  C(R+R  )  , 

c^' 

1  UJ 

(4) 

|o|  < 

|Vtp|  < 

-S|z| 
Ce  '  ' 

1 

'  ^  °  2 

[(e  t  eg 

(5) 

lim 

2  -¥  OO 

\(  j) 

=  X,  X 

»  ( 0 ,  u ) , 

(6) 

lim 

j  >  oo 

=  'f  uniformly 

in 

)(1  -  X/u))’^^ 


on  bounded  subsets  of 


n. 


(PF):  Let  X  .  (n/a,u)  be  given.  For  each  k  >  0  let  (X,ip^)  be  a  solution  of  (PF) 
from  theorem  3.9.  Then  there  exists  an  increasing  sequence  k(j)  +  <»,  j  e  Z* ,  as  above 


and  a  solution  (>,t)  of  (P?  satisfying  (1),  (4),  and  (6)  as  well  as 


. . ,1.2, ' 


(3*)  l9l  -  <  C(«  +  a")“,  where  6 

C^'“ 


(u-xi 


S-m 

4(m-1) 


(Xa-a) 


(B-l) 


Proof!  The  bounds  derived  up  to  this  point  enable  one  to  extract  a  subsequence  of 
solutions  converging  to  the  desired  solitary  wave  solution.  As  In  lenma  4.7  of  Bona, 

Bose,  Turner  (1983),  It  can  be  shown  that  He,  I  Is  bounded  below  by  a  positive  constant 
Independent  of  )c.  It  follows  that  9  >  0.  For  further  details  see  theorem  5.1  In  Bona, 
Bose,  Turner  (1983). 

6.  Example 

In  the  following  specific  example  we  show  how  the  theory  developed  In  the  previous 
sections  Is  applied. 

For  the  primary  flow,  q  ■  (0,V(r)  ,W(r) ),  we  take  W(r)  «  d,  V(r)  «  where  d 

la  a  constant  and  0  Is  an  angular  velocity  parameter,  as  noted  by  Pritchard  (1970)  and 
Benjamin  (1962),  solitary  waves  do  not  exist  when  the  flow  Is  a  rigid  liody  motion, 

V  Or.  In  the  following  calculation  V  -  Or  would  produce  a  linear  equation,  while  we 
require  a  nonlinear  equation.  Since  we  are  looking  for  permanent  waves,  say  of  velocity 
c,  travelling  down  the  cylinder,  we  change  to  a  coordinate  system  moving  with  velocity 
c.  This  renders  the  wave  form  stationary;  W(r)  ”  d  Is  replaced  by  W(r)  »  d  -  c  “  c. 
Using  (2.2), 

r  _  2 

♦  (r)  •  /  sW(s)ds  “ 

0  2  2  2y 

Is  the  stream  function  for  the  unperturbed  flow.  Inverting  this  wo  have  R  ■  — .  Then 

2  ° 

K(r)  -  rV(r)  «  Or  and  In  the  primary  (unperturbed)  flow  l(r)  -  ^  becomes 

2  2  2  ^ 

Kf)  “  Hence  «  120  f  ^  Defining  p(r,c)  by  (2.8)  where 

3  dy  3 

c  c 

\(i(r,z)  •  f(t)  +  yf(r,z),  (6.1) 


lii  Is  the  stream  function  for  the  perturbed  flow,  v  Is  the  perturbed  stream  function 


and  y  Is  a  normalizing  constant  to  be  chosen  shortly.  We  note  that  -  — 

from  (2.9)  we  have 


Y(»  ■  r  ’r  ’zz’  “  [-j"--’*  "  1] 

“  "  R  (IH-y,) 


12(1  .2r2»  c 

[—  •  JTi^ 

(f^-Yy)  [>|I  _  (f+Y.p)] 


-  1] 


12(1^ 


Y9(f+Ys) 


6(1  f  2  ^  2x 
-  -  Y9(r  +  9j* 


0,  so 


Hence 


(6.2) 


We  have  an  equation  ot  the  form 

;(’rr  •  r  *r  ^  ’zz^  * 

6fi^  2  2y  2 

where  h(r«s)  =  0,  X  «  — j—  and  f(r,<p)  »  r  ?  +  ^  9  .  Note  that  X  Is  proportional  to 
2  ®  2 

(1  and  inversely  proportional  to  c*,  the  square  of  the  velocity  of  the  wave  of 

permanent  form  relative  to  the  velocity  of  the  primary  flow.  Choosing  y  by  y  •  c/2, 

we  have  two  cases.  They  correspond  to  c  >  0  and  c  <  0. 

Consider  first  y  >  0.  We  remember  that  in  deriving  (2.9),  (2.10)  we  followed  stream 

surfaces  to  z  -  t  >  and  used  the  knowledge  that  the  Bernoulli  quantity,  H,  and  the 

"circulation**,  I,  were  constant  on  stream  surfaces.  So  the  f  and  h  defined  from 

(2.9),  (2.10)  require  that  t|i  takes  Its  values  In  (0,¥(a)].  In  our  example  this  Is 

(0,ca^/2].  Since  ip(0,x)  -  9(a,z)  «  0,  this  requirement  is  liq>Iied  by 

*  -  9  +  Y9  >0  (6.3) 

*  IT  r 

for  (r,z)  «  (0,a)  x  R.  Using  9^  «  cr  and  y  »  c/2,  we  see  that  (6.3)  Is  equivalent  to 

9^  >  -2r.  (6.4) 


The  restriction  (6.4)  will  certainly  be  satisfied  If  we  can  verify  that  for  our  solution 


(6.5 


Uyl  < 

for  (r,s)  e  Cl  -  [0,a]  x  k.  in  th«  p«rturb«d  flow  the  redial  velocity 

i|i  ye  A 

u(r,c)  - - - - We  know  our  aolutlona  are  poaitlve  (e  “  e)  and  have  a  profile 

2 

In  z  like  e"*  .  In  the  region  to  the  left  of  the  creat,  >  0,  ao  u(r,a)  <  0  and 


we  have  a  wave  of  depreaaion.  Aa  atated  above*  f(r*t) 


r^t  +  t*  -  r^t  +  t^  la,  at 


the  moDent,  only  valid  for  r  c  [0,a]  and  t  >  0  auch  that  TCr)  'f  yt  Ilea  In 


[0,y(a)].  So  (r,t)  muat  aatlafy 


0  <  t  < 


A  2  2 

i  fC*  cr  > 
y  I  2  “  2 


However,  by  extending  f^  ~  t^  aa  an  odd  function  to  t  <  0,  we  aee  that  condition  (H) 

of  aectlon  3  la  clearly  aatlafled  with  Bi“n“2,  -  1,  and  0  “  2/3.  Theae 

c 

condltlona  are  aatlafled  without  any  reatrlctlona  on  the  valuea  of  t  >  0,  ao  our 

extenalon  of  f  ■  r^t  +  t^  outalde  of  0  <  t  <  a^-r^  la  permlaalble. 

Since  the  f  and  h  derived  for  our  exaa^le  aatlafy  (H),  theoraaa  3.7  and  3.9  give 

the  exlatence  of  z-perlodlc  aolutlona  (X,a)  of  aquation  (3.1).  In  (PC)  R  >  0  la 

2  2 

apeclfled  and  In  (PF)  \  «  (0,)i)  ia  apeclfled.  we  know  that  ^  •  x  q,  where  g  «  C 

2 

and  r  <  a,  ao  *  2rg  +  r  g^.  Hence 


<  2r|g|  +  r*(g^|  <  2r(l  +  j)lgl 


From  corollary  4.5  In  Bona,  Boae,  Turner  (1983),  with  N  •  R  In  (PC)  and 

5-m 

f  \i^**”^*  2 

(f  .  -  In  (PF),  Igl  ,  <  C(N+N  ). 

X- 

Hence  (6.6)  will  be  aatlafled  if  (1  a/2)lgl  a  <  1>  Thia  can  be  accompllahed  If  we 

c' 

reatrlct  R  to  R  e  (**>l*g)  ^  ®  or  If  we  reatrlct  X  to  X  e  tX|j,ii)  for 

some  X  €  (0,;i).  With  these  restrictions 

2  60^ 

(i(r,a)  -  T(r)  ♦  Ya(r,E)  -  c(r  +  a)/2,  c  -  , 

has  range  in  I0,T(a)l,  la  a  solution  of  (2.5)  and  gives,  for  each  k  >  0,  a  train  of 


depression  waves.  Since  m  «  2  <  5  and  0  •  a/2a  <  u,  lemma  4.1  of  Bona,  Bose,  Turner 


(1983)  shows  X  <  u  when  R  Is  given.  Inequality  (5.1)  guarantees  that  ip  has 

nontrivial  z-dependence  for  k  sufficiently  large.  Finally,  theorem  5.2  gives  the 

2 

existence  of  solitary  wave  solutions  (X,^)  with  t|»  <•  c(r  t  if)/2.  Since 
2 

|xf^-h^|  =>  Xt  we  know  that  the  deviation  of  the  velocity  fields  of  these  solitary  waves 
from  the  velocity  field  of  the  primary  flow  decays  exponentially  as  |z| 

In  the  case  c  >  0,  the  fluid  Is  moving  downstream  as  viewed  from  the  stationary 
waveform  (W  >  0).  In  considering  the  case  y  -  c/2  <  0,  we  note  that  the  change  In  sign 
of  c  corresponds  to  the  fluid  now  moving  upstream  relative  to  the  waveform  (W  <  0). 
Hence,  the  case  c  <  0  Is  just  the  reflection.  In  z  0,  of  the  case  c  >  0. 

Therefore,  we  again  obtain  a  wave  of  depression. 
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